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PINES METHOD FOR CIRCULAR AND ELLIPTICAL ORBITS 


1 . SUMMARY 

The Pines method for trajectory integration is derived from the equations of 
two-body nwtion. The concept of the perturbation derivative is explained, and 
a formal definition is presented. An algorithm for applying the Pines method 
to numerical integration of trajectories is included, 

2. INTRODUCTION 

The Pines method is a variation-of-parameters method in which initial con- 
ditions on a two-body orbit are treated as elements of the orbit. This 
possible because specification of position and velocity at any time uniquely 
determines a two-body orbit. If at time tg, the position and velocity are 
^ and Vq, at some later time, t, the position and velocity are R and 

V. In the absence of perturbations, the two-body orbit determined by the 
vectors R and V is the same as the two-body orbit determined by the vectors 
^ and If perturbations are present, however, R^ and must 

gradually change so that the two-body trajectory on which they lie will inter- 
sect the position R with velocity ^ and are the elements of 
the osculating orbit, and the rate at which they change is called the 
"perturbation derivative" of R^ and V^. If the perturbation derivatives 

of R„ and are known, R^ and can be calculated for any time by 
— u ^ ~u 

numerical integration. Once ^ and ^ for a given time are known, the 

“HJ 

instantaneous two-body orbit is known and the actual position and velocity 
(R, V) can be calculated from the formula* for two-body motion. This is 
the essence of the Pines method. 

3. THEORY OF TWO-BODY MOTION 

The formulas for two-body motion are well known and may be found in a number 
of references, including reference 1 (chs. 1 and 2) and various astronomy texts 


(such as ref. 2 (ch. V, p. 107)). For the convenience of the reader, however, 
derivations of the formulas introduced in this section have been included in 
the appendix. 

The position and velocity vectors for some point on a two-body orbit are given 
by eqs. (96) in the appendix 


R = 1 - ~ (1 - cos q) ^ 


( 1 ) 


+ /t - t - ^ V 

V “ ^ r 


^ sin q ^ + 1 - -^1 - cos q) 


( 2 ) 


where quantities with subscript zero are evaluated at lime tg, q is the 
difference in eccentric anomaly (E - Eq), m is the gravitational constant 
times the mass of the central body, a is the semimajor axis, and r is the 
distance from the central body. Giver, an initial position and velocity, we 
may find position and velocity at some other time by using equations (1) and 
(2). We obtain a by using eq. (101) in the appendix 



(3) 


where 


V 


0 



(A\ 

\ • s 


2 



and 



( 5 ) 


The quantity q = E - Eq may be obtainad by solving Kepler's equation 
(eq. (93) in the appendix) in the form 


Eqs. (1) and (2) are valid whether we go forward or backward in time from the 
initial point. We could express position and velocity at the earlier time in 
terms of position and velocity at time t and let this earlier time be tg. 
The result is identical to eqs. (1) and (2) except that subscripted and 
unsubscripted quantities are interchanged and q is replaced by -q. 

Before we do this, let us introduce the scalar functions f and g given by 


f . 1 . Ml . cos q) 
^0 


and 



Differentiating eq. (7), we find 


f - a „ dE 


(7) 


( 8 ) 


(9) 
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We are confining ourselves to two-body motion for the nKwnent; therefore, a - 
and rQ are constants. From eq. (84) In the appendix, we have 


So eq. (9) may be written as 

f . .-JE sin q (11) 

"o 

Differentiating eq. (8) yields 

9 - ’ 

And by applying eq. (10), we find that 

g = 1 - - cos q) (12) 

Using eqs. (11) and (12), we may rewrite eqs. (1) and (2) In the form 

R-fRo + gV(,l 
1 - ^Bo"9Voj 

where, after collecting our formulas together, 

f = 1 - ^1 - cos q) (from 
*^0 


(13) 


we have 
eq. (7)) 
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(from eq. (8)) 


. t - t. - SLz^ 



f * - sin q 
0 . 


(from eq. (11)) 


g = 1 - p (1 - cos q) 


(from eq. (12)) 


If we invert eqs. (1) and (2) to express ^ and in terms of R and V, 
we have 


^ = [l _ i{l _ cos q)] R - 


1 ^ ' 

(t . t„l i 

I - 

“ m 


V = 

■H) rr. 


sin q_R + 1 - ^1 - cos q) V 
*^0 


or, comparing with eqs. (7), (8), (11), and (12), we find that 


^ = gR - g V 


^ = -fR + fV 


(14) 


If we rewrite Kepler's equation (eq. (6)) in the same way, we have 



(1 - cos q) - 



q 


(15) 
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Note that Instead of ^* 1 q * ^ now have ^*V; and Instead of fg, we have 
r. We can use either eq. (6) or eq. (15) to eliminate the t1me> dependent 
term In g. If we use eq. (6)» we have 


9 = y - cos q) + sin qj (16) 

and If we use eq. (15). we have 

9 ■ “ - cos q) - r^aiiT sin ql (17) 


*;1ther expression for g may be used - whichever is the most convenient. 
From equation of the orbit (eq. (83) In the appendix). 


r = a(l - e cos E) 

= a[l - e cos(e - Eg + Eg)] 

. a(l - e cos Eg cos q + e sin Eg sin q) (18) 

From the appendix (eqs. (90) and (92)), 


e cos Eg . 1 - 

(19) 

e sin Eg • ^ Rg-Vgy>/:r 

(20) 


Substituting into eq. (18) yields 

r = a(l - cos q) + rg cos q + Rg*V^ ^ sin q (21 ) 
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Eq. (21) allows us to calculate r If we know only the Initial conditions and 

q- 


In the next section, we will need an expression for R-V In tenns of the 
initial conditions. This Is found by taking eq, (90) from the appendix In 
the form 


R*V. *^ua e sin E 


and using the following expansion*. 


e sin E ® e sin 


(' - ^0 * ^o) 


e cos Eq sin q + e sin Eq cos q 


(.•?)■» 


11 - jS) sin q + cos q 
* ' ^ 


where we have used eqs. (19) and (20). With this result, we may write 


R.V =^ua ^1 - sin q + Rq*^ cos q 


( 22 ) 


The formulas presented in this section allow us to obtain the position and 
velocity of an object moving along an elliptical orbit provided that we have 
the initial conditions. We would proceed as follows: From the initial con- 
ditions ^ and Vq, we can calcualte rQ, Vq, and Rq*V^j- Knowing rQ 

and Vq , we can calculate a by using eq. (3). We can then solve 
Kepler's equation in the form of eq. (6) to find the value of q. Eq. (21) 
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can be used to find r, and we can then substitute the values of a, 
r» and q Into the formulas for f and g and evaluate R and 
V by using eqs. (13). 


4. THE PERTURBATION DERIVATIVE 

In the presence of perturbations, the vectors ^ and V^ will change slowly 
in time. In order to calculate the rate of this change, we shall define a 
“perturbation derivative," which Is the time rate of change of a quantity with 
ivspect to what It would have been In the absence of perturbations. 

Let Q be some function of time along a real trajectory. If we expand Q 
It Taylor series about the time t, we have 


2 

q(t + At) • Q(t) + 3^t + 7 ^At)^ .... 

dt 

If the rate of change of Q along a two-body orbit at time t Is written as 
Q^, then the perturbation derivative Is Qf * ^ 



■QjtQ 


T 


The T subscript specifies the perturbation derivative, and the t subscript 
specifies the derivative for the Instantaneous two-body orbit with the same 
state vector. 

To first order, then, we may write 

Q(t + At) - Q(t) + Q^At + Q.|^t (23) 
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The sum of the first two terms on the right Is Just the value that Q would 
have at t + At If It were evaluated along the two-body orbit. We abbreviate 
this by 

Q*(t + At) - Q(t) ♦ Q^At (24) 

Solving eq. (23) for Qj and taking the limit as At approaches zero, we 
have 



Eq. (25) will be our formal definition of the perturbation derivative. 

The perturbation derivatives of R and V can readily be calculated from 
our definition and the equation of motion for two-body motion, which is 

, uR 

( 26 ) 

r 

Expanding R In Taylor series yields 

R(t + At) = R(t) + V(t)At + 7 V(t)(At)^ ■»■... (27"^ 

V(t) is the instantaneous two-body velocity; therefore, 

R(t) + V(t)At » R*(t + At) (28) 
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SmODTCTOV i >Y 

To first order, 

R(t + At) ■ R*(t + At) (29) 

Therefore, applying eq. (25), we find that 

- 0 (30) 

The equation of motion, eq. (26), Is just the statement that 

dV uR 

^ - -g-+ F (31) 

r 

The first term on the right Is the two-body acceleration. Expanding 
V(t + At) about t yields, to first order. 

uR 

V(t + At) « V(t) - At + FAt (3?) 

- 

A 

but the first two terms on the right are just V (t + At). We have 

V(t + At) « V*(t + At) + FAt (33) 

Therefore, by our formal definition. 
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We have now obtained the Important results 


(34) 


Rj ■ 0 (from eq. (30)) 

and 

■ £ (from eq. (34)) 


These are valuable results because any quantity for which we might need to 
calculate the perturbation derivative can be expressed as a function of R 
arc' V. Taking the perturbation derivative of eqs. (14) and using eqs. (30) 
{^nd (34) results In 


(Bolj ' SjE •• 9# - 9l 
(»o)t ■ 


(35) 


Expressions vor f^. g^* fT. and g^ are needed; but before they can be 
derived. It Is necessary to resolve a certain ambiguity about what Is meant by 
Initial conditions. 


One way of defining our Initial conditions would be to say that tQ Is a 
fixed time and that and Vq are the position and velocity at tp for a 
two-body orbit determined by R and V at time t. Under this assumption, 
t - tQ Is just the elapsed time for a two-body orbit and 

(t - • 0 (36) 
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This choice leads to the set of equations originally presented by Pines 
(ref. 3). 

Another way of defining initial conditions would be in terms of the eccentric 
anomaly. If Eq is fixed and ^ and are the position and velocity on 
a two-body orbit determined by and ^ when the eccentric an^aly is E, 
then q = E - Eq is just the change in eccentric anomaly along a two-body 
orbit and 

,T = (e - = 0 (37) 

This choice yields the set of equations presented by Pines in ref. 4 and is 
the choice used in the following work. 

5. EQUATIONS OF MOTION 

To save writing, let 


<‘o = ^0-^ 


d = R-V 


Eq. (22) may be written as 


d sin q|l - ^ 


(38) 


since d = R»V, d^ = or by using eqs. (30) and (34), 


d 


T 


» R«F 


(39) 
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The perturbation derivative of a can be calculated by using 


1 1 



(40) 


Taking the perturbation derivative of eq. (40) and using eqs, (30) and (34) 
yields 


aj = 2a‘V«F/y 


(41) 


In a similar manner, all the perturbation derivatives of f, g, f, and g 
can be calculated by using the assumption q^ = 0. Applying the perturbation 
derivative to eq. (7), we have 


f j = I a(l - cos q) 


Iro): 


/r. 


(42) 


Taking the perturbation derivative of eq. (17) yields 


9t- 


r sin q d(l - cos q) . a(l - cos q) j 


(43) 


Direct application of the perturbation derivative to eqs. (11) and (12) gives 


• • 


"T 

I? 


Oil 

r 


0 J 


(44) 
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f 


and 


9t 


_ _ 1 > cos q 


(45) 


We now have expressions for everything in eqs. (42) through (45) except 
(rgly. If we invert eq. (21), we have 


rQ = a(l - cos q) + r cos q - sin q 
and direct application of the perturbation derivative yields 

fo)r " " Vf (^T ^ 9 (46) 

where we have used the fact that 

>•7 ' [(i’i) JT " " ° 

We now have everything we need to calculate the derivatives of ^ and \L 
for any value of q. We could obtain q by solving Kepler's equation; 
however, it is convenient to introduce a seventh variable of integration from 
which q may be calculated. 

Let W = q VT The total derivative of W is 


dW 

3t 



2 ^ 


da 
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From eq. (10), 





Eq. (10) applies because the perturbation derivative of q is zero. The 
rate of change of q is the same as it would be for a two-body orbit; there- 
fore. we can use eq. (10) to express it. 

Since a is constant for a two-body orbit, the total derivative of a equals 
the perturbation derivative ^ = a.j. and we have, using eq. (41), 


^ + ^ VF 
dt r y - - 


(47) 


We can integrate together with and «t each step and 

calculate q by using 


q = W/yr (48) 

Collecting together the differential equations which are to be integrated, we 
obtain the following set of equations of motion: 


^R^.gTR-grV-gF 


(from eq. (35)) 


^ ^ ^ VF 

dt r y -- 


( 47 ) 


where we have used the fact that the total derivative of the initial con- 
ditions is equal to the perturbation derivative. 

We now have everything we need to integrate along the trajectory. 
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6. INTEGRATION ALGORITHM 


STEP 1: Initialize variables 

At the beginning of the Integration, set 



STEP 2 : Function evaluation 

The following order of evaluation of the variables may be used to calculate 
the derivatives of the Initial conditions and the auxiliary variable W. 



q = W/Va“ 




f = 1 - a(l - cos q)/rp 
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r » aO - cos q) + Tq cos q + /Tsinq 
f = - sin q 
g “ 1 - f 0 - cos q) 
i = fBo * 9Yo 

V = fBo + 9J!o 

F “ £(R, V) (Evaluate perturbations.) 


dj = i*£ 


d =Vvia y - sin q + d^ cos q 

(’•q)t ' (’ - 9)aT -Jf * 5t) 


f_ - ^(1 - cos q) 




^ Bfl = jji - 9jl - gL 

^ = JC+ *M v-f 
at r y — 

Note that the total derivatives of ^ and are just the perturbation 
derivatives because they are constants along the two-body orbit, whereas W 
does change along a two-body orbit, it being a function of eccentric anomaly. 

STEP 3: Integration 

After having obtained the needed derivatives, new approximations of R^, V^, 

and W may be obtained according to the integration scheme being used. We 
then return to step 2 for the next function evaluation and iterate. 

7. USE OF INITIAL TIME AS AN ELEMENT 

The form of the equations of motion as presented in section 5 have one 
peculiarity, the fact that all the quantities being integrated are slowly 
varying except the seventh variable W which varies rapidly along a conic. 

The variable W was introduced to provide a way of calculating the difference 
in eccentric anomaly and to avoid having to solve Kepler's equation. For 
near-circular orbits, however, the solution of Kepler's equation is not time 
consuming, therefore the use of an alternate form of the Pines method should 
be considered. 

The value of q may be found by solving eq. (6). The left-hand side of 
eq. (6) is the mean motion 


M 



(48) 
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and eq. (6) rnay be rewritten In the form 


M - q + i (1 - cos q) - (1 - sin q (49) 

^ * 

A Newton-Raphson technique may be used to solve eq. (49). The Initial 
estimate of q is 


% ' M 

and successive approximations are found according to the scheme 


(50) 


M - 
d?”n 


where 



+ (1 - cos q ) - (1 

Npa 


r> »i" % 


(52) 


and 


d 


M 


n 




(53) 


For elliptical orbits, the eccentricity lies between 0 and 1. The Newton- 
Raphson solution of eq. (4) has been found to converge within five iterations 
for eccentricities between 0 and 0.5. Once q has been found, the Pines 
method proceeds as previously outlined, except that the differential equation 
for W is removed and replaced by a differential equation for tg. 
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Tha Initial time tQ Is not a constant because of our definition of Initial 
conditions. Because our Initial conditions were defined In terms of eccentric 
anomaly so that eq. (37) would be satisfied, eq. (36) Is not satisfied, there- 
fore tjj Is not fixed, but varies slowly as the other Initial conditions. 

A differential equation for may be most easily found from eq. (15) 
which can be rewritten as 

t - to q - ^ <10 - cos q) - j) sin q (54) 

The total derivative of t - tp is obviously 

^ " Hr 


but it is also true that 


- tg) = (t - ^ (t - 1„), (56) 

where, as before, the t subscript means the time rate of change of the 
quantity when initial conditions (in this case, tg) are held constant and the 
T subscript indicates the perturbation derivative. 

Now 

(t - = 1 

therefore eq. (56) is equivalent to 
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(57) 


- tj|) • 1 + (t - 

Using eqs. (57) and (55), we find that 

^ - - (t - tp)^ (58) 

which means that the time rate of change of tQ Is the negative of the 
perturbation derivative of the right-hand side of eq. (54). 




3(q - sin q) + r sin q/a 


- |r 0 - cos q)]a, - -J” (^ - cos q) d ' 

Vwa j 


(59) 


where and are the perturbation derivatives of a and d (given by 
eqs. (41) and (39), respectively). Eq. (59) replaces eq. (47) and the seventh 
variable is now t^^ Instead of W. 

If the Pines method is set up in the form described in this section, the 
technique will provide an exact solution to the two-body problem. When no 
perturbations are present, all perturbation derivatives are zero and a state 
vector may be propagated ahead to any time to the limit of machine accuracy in 
a single step. 


8. CONCLUSIONS 

The Pines method for circular and elliptical orbits outlined in the previous 
sections provides a more stable and accurate method for integrating the equations 
of motion than does the Cowell method. The numerical accuracy of this method 
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Is not as great as that of the KS method Cref. 5), but this method Is easier 
to program, requires less storage, and has a shorter cycle time. Only seven 
simultaneous differential equations must be Integrated with the Pines method, 
whereas the KS method requires 10. Time Is the Independent variable In the 
Pines method, whereas time Is an element In the KS method; thus, the Pines 
method Is considerably simpler to program In typical applications. Test runs 
with this method using an eighth-order potential and Including perturbations 
due to the sun, the moon, and drag have shown the error to be 11.5 meters 
after 2 days when a time step of 1/50 revolution for a Shuttle-type orbit Is 
used. 
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APPENDIX 


TWO-BODY MOTION 

CONSERVATION OF AN6UUR MOMENTUM 

The equations of motion for a point mass In an Inverse square attractive 
force field are, In vector form, 

R + Vi*0 (60) 

r’’ 

Forming the cross-product of R with eq. (60), we have 

R X R - 0 

which is equivalent to 

^(R X « 0 

from which we conclude that 

R X ft ■ h (61) 


where is a constant vector. It can be seen, by forming the dot product 
of R with eq. (61), that R*h * 0; consequently, the motion of the orbiting 
body is confined to a plane that Is perpendicular to h and that contains 
the center of attraction. Since this is true, it is convenient to choose a 
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coordinate system In which the x and y axes are In the plane of motion 
and In which the z*ax1s Is parallel to h. Writing out the equations of 
motion In scalar form yields 


X X - 0 

r 

Let r, e be polar coordinates such that 


( 62 ) 


X ■ r cos e 
y » r cos e 

Differentiating twice and substituting In eq. (62) yields 


(63) 


r cos 6 - re sin 6 - re^ cos 6 - 2re sin 0 ♦ cos 9 - 0 (64a) 

r 

and 

r sin 0 - re cos 9 - r6^ sin 6 + 2r9 cos 9 + ^ sin 0 « 0 (64b) 

r 

Multiplying eq. (64a) by cos 0, multiplying eq. (64b) by sin e» and adding 
yields 

r - ^ = 0 (65a) 

r 


Multiplying eq. (64a) by sin 0, multiplying eq. (64b) by cos 9, and sub- 
tracting yields 


24 



-re - 2re ■ 0 


or 


^(r^e) • 0 (65b) 

Eq. (65b) implies that re is a constant; and if we compare with cq. (61) 
written in polar coordinates, we see that the constant is the magnitude of ji. 
Therefore, 

r^e * h (66) 

Using eqs. (65a) and (66), we find the equation of the orbit by making the 
change of variable r = 1/u and using eq. (66) to eliminate the time deriva- 
tive. Set r = 1/u. Then, differentiating with respect to time yields 



1 du 
'T d0 


From eq. (66), however, 


• 0 

e = h/r 



Therefore, 


du 

de 


Differentiating with respect to time again yields 


k ^^u X 2 d^u 

r = -h — y 0 " -h u — T 

de^ 


(67) 


( 68 ) 


( 69 ) 
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Substituting for r and 0 In eq. (65a) yields 


d u p 

—K + u - -j- 

de'^ h 


( 70 ) 


The solution of this differential equation Is easily found and can be written 
in the form 


u = ^ j^l + e cos(0 - w)j (71) 

where e and w are the constants of integration. Replacing u by 1/r 
yields the equation of the orbit in polar coordinates: 


' T + e co?(0 - 0 )) (72) 

The constant e Is called the eccentricity, and w is the longitude (or 

p 

argument) of periapsis. When 0 = w, r has a minimum value. Let P = h /y, 
where P is called the parameter (or semilatu s rectum ). If the coordinate 
axes are rotated so that the x-axis is along the line of periapsis (line 
apsides) and we define v = 0 - w, where v is called the true anomaly, then 


X = r cos V 
y = r sin V 


Eq. (72) can be written as 


( 73 ) 


r + er cos v = P 
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or 




X + y + ex = P 




P - ex 


(74) 


Squaring eq. (74), we have 


2 2 2 2 2 
X + y = P - 2ePx + e x 


(75) 


If we collect the terms containing x together and complete the square, we 
have 



which can be put in the form 


where 



(76) 


(77) 


(78) 
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provided that e < 1. If e < 1, eq. (77) Is the equation of an ellipse with 
the right-hand focus at the origin of the coordinate system and the other 
focus on the -x axis. The geometric center of the ellipse will be at the 
point (-C, 0), as can be seen In figure 1. If we draw a circle of radius 
a about the center of the ellipse and then draw a line through the point 
(x, y) parallel to the y-axis to Intersect theicircle at point P, the angle 
subtended by the arc of the circle from point P to where the circle cuts 
the x-axIs Is called the eccentric anomaly, E. By Inspection of figure 1, 
we see that 

X = a cos E - c = a(cos E - e) (79) 

where we have used the fact that c = ae from eqs. (78). By substit uting 
with eq. (79) Into eq. (76) and replacing c by ae and b by P/‘^lT~-i?, 
we obtain 


y = . sin E = Sin E (80) 

yJTJ 

So, we may write 

^ A 

R = a(cos E - e)x +\aP sin E y (81a) 

and, by differentiating, 

V = -aE sin E x E cos E y (81b) 


^ A 

where x and y are the unit vectors along the x and y axes. 
Using eqs. (81a), (81b), and (61), we find that 

h = hz = (a^P)^(l - e cos E)Ez 
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Therefore, 


(a^P)^(l - e cos E) 


But 




Figure 1.- Elements of the orbit. 
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So 




a(l - e cos E) 

Now r * so If we use eqs. (79), (80), and (78), we find that 


(82) 


r = a(l - e cos E) 


(83) 


Therefore, 


E = l# 

r^a 


Substituting back into eqs. (81a) and (81b), we have 


A A 

R = a(cos E - e)x +“yaP sin E y 


V. = “ f ^ap sin E X + cos E y 


(84) 


(85) 


By solving eq. (85) straightforwardly with the use of Kramer's rule, we 
obtain 


X = 


cos E 


R sin 

- Xu 


E V 


( 86 ) 


y "MF 


a sin E ^ ^ ^ 


r - 


■y/J 


KEPLER’S LAW 


We need an equation for E, the eccentric anomaly. Using eqs. (84) and (83), 
we find that 


30 


( 87 ) 


(1 - e COS E) ^ 

which is easily integrated and yields 

^(t - tg) - E - Ep - e ^sln E - sin E^j (88) 


which is one form of Kepler's equation, where E is the eccentric anomaly at 

time t and Eq is the eccentric anomaly at t^. 

We would like to get eqs. (85) and (88) into a form such that eccentric 

anomaly only appears in the form E - Eq. We can write 


sin E = sin^E - Eq + Eq| 

fn(E - Epj + sin Ep cos(e - Ep) 


= COS Eq si 


and put eq. (88) in the form 


- ‘o) “ ^ - ^0 ^ ® ^o[' ■ - ^o)] ■ " ^0 - ^o) 

2 

We can use eq. (85) to calculate If we replace P by a(l - e ) 

(eqs. (78)) and use eq. (83), we find that 

R*V 

= e sin E (90) 

Evaluating eq. (90) at time t = t^ gives 



e sin E, 


(91) 
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Evaluating eq. (83) at t« and solving for e cos E yields 

^ 0 

(| - r)' ® ^0 

Eqs. (91) and (92) can be used to put Kepler's law In the desired form by 
substitution Into eq. (89). 



POSITION AND VELOCITY FROM 
INITIAL POSITION AND VELOCITY VECTORS 

Eqs. (86) expressed the unit vectors parallel to the semimajor and semiminor 
axes in terms of position and velocity vectors. Since the unit vectors are 
constant vectors, the relations (eqs. (86)) hold for any time, t^. 



Direct substitution with eqs. (94) into eqs. (85) yields 
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(95) 


- ‘I ■ ^o) ■ ® 

- Eo)+ e s’" Efl - ® ^]{j!o 

- “ ■ ^ s1n(E - E JEo + |f [cos (e - Eg)- e cos 



Eqs. (95) have two undesirable features; They contain terms in eccentric 
anomaly that do not depend only on E - Eq, and the eccentricity appears. 

Using eqs. (83) and (92), we can get rid of the terms e cos Eq and e cos E. 
From Kepler's law (eq. (88)), we can write that 



and use this expression to eliminate the terms e sin E and e sin Eg. The 
result is 



Eqs. (96) are the usual form for the f and 'g expansion for two-body 
motion. 


33 



ORIGW^L P^Gb 1 ■ 


CONSERVATION OF ENERGY 

Up to this point, we have everything we need to describe two-body motion 
except a formula for a. In order to obtain such a formula, we need to derive 
the conservatlon-of-energy formula. This Is done by dotting equation (60) with 
V and using R ■ 


r'^ 

Noting that r * (R*R)*®, we can put eq. (97) In the form 

rt(¥)- ^ ■ 0 

Therefore, 


5^ - ^ = K (98) 

where K is some constant. We can determine the value of K by evaluating 
eq. (98) at any point on the orbit. Let E * 0 and use eqs. (85) to get v 
and r. We find that 



Tq = a(l - e) 

And substituting Into eq. (98) yields 


uP 

2^V~^~e7 



K 


(99) 
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But P 



(eqs. (78) )t so we have 


V 

a 



2(1 - e)^ 



K 


which reduces to K 



Eq. (98) can now be written as 



( 100 ) 


which Is the familiar Vis Viva Integral. From eq. (100), an expression for 
a Is obtained. 

»=l/(|-J^) (101) 

COLLECTION OF FORMULAS 

The results from this appendix that we shall need are as follows: eqs. (96), 

R = jl - ijl - cos(e - Eq)]Ir, 

V = .^s1„(e.Eo)Eo- jl -f[l-cos(E-Eo)]jyg 
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Kepler's Law in the form of eq. (93), 


R.-V, 


- ‘o) = (e - Eo) * ^ [l - COS 


- 1 


(■ ■ ^o) 


and the relations 



(eq. (90)) 

(eq. (92)) 

(eq. (84)) 
(eq. (101)) 


r =» a(l - e cos E) (eq. (83)) 
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